function run
lambda = 2*pi^2;
ep = 0;
nx = 30;
ny = nx;
tol = 10e-5;

%gen guess
guess = ones((nx+1)*(ny+1),1);

    [jaco]=gen_j(guess, lambda, nx, ny);
    
    [R] = gen_ar(guess, lambda, ep, nx, ny);
    sol = jaco\R;
% calc
%[sol, it] = newton(guess, lambda, ep, nx, ny, tol)

%plot

plot_sol(sol, nx,ny);
det(jaco)